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We study semi-classical slow light propagation in trapped two level atomic quantum gases. The 
temperature dependent behaviors of both group velocity and transmissions are compared for low 
temperature Bose, Fermi, and Boltzman gases within the local density approximation for their 
spatial density profile. 
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Recently, a dramatic demonstration of slow light prop- 
agation, down to 17 (m/s) was reported 0,0. Achieve- 
ment of such an extremely low speed owes mostly to the 
widely discussed phenomenon of electromagnetically in- 
duced transparency (EIT), which makes propagation of 
light in an otherwise opaque medium possible ||. In 
related theoretical and experimental studies, three level 
atomic vapours are the typical medium. Although Bose- 
Condensation is not crucial to EIT or slow light propa- 
gation [jl]||] , the long coherence time of such a degenerate 
quantum medium does prove to be advantageous. Simi- 
lar observations have since been carried out in optically 
dense hot rubidium vapours [~ 90 (m/s)] [|J and in a 
optically thick Pr:YSO crystal [~ 45 (m/s)] f§. Possi- 
bilities of vanishing or even negative group velocities, in 
a coherently driven Doppler broadened atomic vapours 
have also been discussed 

Ultra-slow light propagation offers many potential 
applications since macroscopically the medium can be 
viewed as possessing a high index of refraction, albeit 
within a narrow spectrum. The accompanied super- 
high nonlinear coupling between weak fields in a long 
coherence medium such as a Bose-Einstein condensate 
(BEC) opens up novel regimes of quantum nonlinear op- 
tics M. Application to quantum networks and quan- 
tum information processing, including quantum entan- 
glement of slow photons || , non-classical (e.g. squeezed) 
and entangled atomic ensembles quantum memo- 
ries pj, have been proposed. Their implications to 
quantum non-demolishing measurements and high pre- 
cision spectroscopy using squeezed light have been sug- 
gested through enhanced acousto-optical effects 111 ] and 
narrow-band sources for non-classical radiation |12 . 

The aim of this paper is two fold: (1) Inspired by 
the recent theoretical modeling of slow light propaga- 
tion in BEC by Morigi and Agarwal we ask 
the question of comparative differences related to spa- 
tial density profiles of different quantum gases; (2) We 
explore a simpler model composed of two level atoms 
f[4} |. In our formulation, we consider the semi-classical 
propagation of a laser pulse through an ultra-cold quan- 
tum gas of two level atoms described by a density profile 
p(r, t). At very cold temperatures atoms are highly delo- 
calized and continuum treatment of the g medium 
is desirable. The polarization density operator is thus 



P(r,t) = ct ge i/)l{r,t)il>e(r,t) + h.c. @ with d ge the elec- 
tronic dipole transition moment between ground state 
\g) and excited state |e). Tjj ri=gte (r,t) is the second quan- 
tized field operator for atoms in state \g) and |e) respec- 
tively. An incident laser pulse polarizes the medium and 
the induced dipoles of atoms subsequently emit coher- 
ent radiation. For most of aforementioned applications, 
the atomic density profile remains an essential ingredi- 
ent. In addition, careful determination of the optical 
response inside a dense medium requires the inclusion 
of local field corrections. In case of two level atoms as 
studied here, it is necessary to also have a significant 
transmission through the medium. At any point within 
the medium, the electric field is given by the superposi- 
tion of the incident field and the radiated field, the latter 
excludes the self-field radiated by the dipole at this same 
point. Adopting the standard approach, we write the to- 
tal local field acting on the atoms Ep = Ep + Ep + Ey , 
as a sum of the incident field (Ep), the field generated by 
other atomic dipoles (Ep) and the field due to vacuum 
fluctuations (Ev)- The dipole field depends on the type 
of medium as well as its density profile. 

Following standard procedures we eliminate fast 
radiation field dynamics. It is then possible to recognize 
that Ep is given by [Pt], the retarded transverse polar- 
ization field We obtain 
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where V/ {r\ indicates that integral over the interaction 
region excludes a small region around f. Proper evalua- 
tion of the integral in the vicinity of r yields the Lorentz- 
Lorenz correction (local field correction) p7| , 

An - 

E L = E+—P. (1) 

The field E, including contributions of all dipoles in the 
interaction volume V, is the macroscopic field governed 
by the Maxwell equations. It is important to point out 
that the Lorentz-Lorenz relation Eq. (|l|) holds true for 
both linear and nonlinear media ||18|] . To maintain suffi- 
cient transmission, we assume the incident field (of cen- 
tral frequency u>) to be far-off resonant (from atomic res- 
onance ujq) with a large detuning A = lu — loq. When 
the incident pulse is weak, the linear response theory 
probes the unperturbed density profile p(r, t) sa p(f) = 
4>l(r, 0)ip g (r, 0). For stronger higher input fields, the den- 
sity of the gas needs to be determined self-consistently 
Jl9[ . The linear optical response function (the electric 
susceptibility, \) follows the Clausius-Mossotti relation 



ap(r) 
'-ap(f) 

involving Lorentz-Lorenz shift as a first correction 
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5|,|19[-|24|. a = | dge | /h A is the approximate off- 



resonant atomic polarizability and 7 is the excited state 



spontaneous emission rate. Higher order quantum cor- 
rections to this formula arise from multiple scattering of 
photons within neighboring atomic pairs as well as from 
many-body correlations. These effects are neglected in 
this study for the coherent propagation effect. Recent 
studies of these corrections found the weak field optical 
response Eq. (Q) works well in available experimental 
density regimes for both Bose and Fermi gases 

p4| , p5| . The main physical reason allowing for such a 
simplification is the experimental effort in assuring a low 
optical density x(^*) = ap(r) to provide reliable probes of 
the system. Local field correction is significant when the 
number of atoms in a characteristic volume 

is comparable to unity. Here, fez, is the wave number of 
the near-resonant incident field with detuning a A (« IO7 
in this study) from the transition wavelength A = 589 
(nm). This gives V a « 4 x 10~ 15 cm 3 . Hence, only for 
densities around f 14 cm -3 or higher, will local field start 
to cause qualitative/quantitative modifications to the op- 
tical response of an atomic cloud. For a non-interacting 
condensate of 10 6 atoms, it is estimated [ p"3"| using trap 
parameters of Ref. |Q that an equivalent homogeneous 
density is p w 8 x 10 15 cm -3 . The presence of a re- 
pulsive interaction between atoms lowers that value to 
p « 3 x 10 13 cm -3 (using Thomas-Fermi approximation). 
Therefore, it is more serious to neglect local field effects 
in a BEC model of non-interacting atoms than our semi- 
classical treatment for an interacting condensate. Al- 
though the above estimates could justify the neglect of 
Lorentz-Lorenz shift, we decided to keep it for a couple 
of reasons. First, these estimates rely on a homogeneous 
density profile. Yet, the density of an actual condensate 
can be much higher near the trap center. For an inhomo- 
geneous condensate the group velocity is determined by 
a spatial averaged optical response, it is therefore neces- 
sary to check consistently any deviations from the column 
density approach used for a homogeneous density distri- 
bution. Second, considerable arguments (see discussions 
in next paragraph) point to the fact that local field ef- 
fects, many body correlations, and multiple scattering 
effects are of the same order. Keeping Lorentz-Lorenz 
shift thus allows for a consistent check against approxi- 
mations in neglecting the many body (e.g. second order 
density-density) correlations as well as multiple scatter- 
ing. 

It is a challenging task to include quantum many-body 
correlations and multiple scattering of the medium into 
an effective optical response (effective refractive index). 
In an earlier attempt, a perturbative density expansion, 
similar to the Virial expansion, was employed in Ref. 
pof to calculate the first two corrections. The effect was 
found to be negligible for atomic densities <C k\. Later, 
it was shown that such an approximation is equivalent 
to truncating an infinite hierarchy of equations for cor- 



relation functions 21 22|,|25j. Furthermore, it was shown 



under the truncation approximation that correlation cor- 
rection terms are of the same order as the local field cor- 
rections. A notable exception is the low density Fermi 
gas in a Cooper paired BCS state pq |, a case not consid- 
ered here. The neglect of quantum correlations even in 
the ideal fermi gas case as studied here may affect some of 
our results. Other recent studies include self-consistent 
analysis in [ p^|jis|j2^| and a Dyson equation formulation 
in g]. 

The real part of the susceptibility (%') is related to the 
index of refraction from which we can calculate the group 
velocity according to j| 

"9( u . r >— J- ~ SjF» ( 4 ) 

where c is the speed of light in vacuum, and the probe 
pulse central wave number is k — uj/c. In obtaining the 
above equation, we have neglected a small contribution 
from dx'/dk, and assumed the susceptibility changes 
slowly over an optical wavelength 2n/k. The result is, 
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In the reported measurements of slow group velocity 
two separate parameters are needed: 1) the effective 
propagation length of the medium; 2) the delay time of 
the probe pulse. The ratio of these two quantities defines 
the effective group velocity |J. We adopted the defini- 
tions as used by Morigi and Agwaral jl3| for a three level 
Bose gas, and appropriately generalized to the case of an 
interacting Bose gas. We then compare with correspond- 
ing results for an (ideal) Fermi and Boltzman gas. 

The ground state atomic density profile for a Fermi 
gas is obtained following the semi-classical approximation 
given by 

where (3 = l/(fcsT). The Hamiltonian H(p,f) describes 
ground state atomic motion inside an external trapping 
potential V(r) = Mu 2 (r 2 + e 2 z 2 )/2. Lu r is the radial trap 
frequency, e is the aspect ratio, and M is the atomic mass. 
The probe light is assumed to propagate along the (long) 
z— axis of the atomic cloud. The chemical potential fi is 
determined by the normalization N = J drp F (f) with TV 
the total number of atoms. In the low temperature limit, 
p can be described within the Sommerfeld expansion, 
while the high temperature behavior is found by direct 
integration similar to a Boltzmann gas p7| , 

_ / - In [6(T/7» 3 ]//3, T > 0.55T F , 
M ~ \ E F [l - ir 2 (T/T F ) 2 /3], T < 0.55T F , {<) 

where Tp and E F are the Fermi temperature and Fermi 
energy. 

For an interacting Bose gas, the ground state density 
profile is computed using the analytic fitting function 



as developed in |j28| . Below the condensation temper- 
ature Tc, the condensate component is described by the 
Thomas- Fermi approximation [Eq,E9|, 
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with T c is the critical temperature. U — ATTh 2 a sc /M with 
a sc the atomic scattering length. 0(.) is the Heaviside step 
function, g n (x) = x J /j n , and At is the thermal de 
Broglie wavelength. The chemical potential at high tem- 
peratures is determined by solving Lis(z) = (T/T C )~ 3 C(3) 
in terms of the fugacity z = e^ M . Lia(.) and C(3) are the 
third order polylogarithmic and Riemann-Zeta functions 
respectively. At low temperatures, it is found that |28|, 
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with /iTF the Thomas-Fermi approximation chemical po- 
tential. The condensate fraction is then given by 



Ao _ /TV _ C(2) (T_ 
N \T C ) V ((3) \T C 
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with a scaling parameter r\ defined as 

KB-Lc 2 V a hoJ 
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(10) 



(11) 



a\ w = ^Jfi/ (M e 1 / 3 ^,.) is the average harmonic oscillator 
length scale p0| . 

Finally, for a classical gas we use 



Pc(r) 



(12) 



with the chemical potential /i = - In [6(T/T F ) 3 }/ /3. 

From each of the above three distributions for p(r), we 
first determine the effective length of the medium using 
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(13) 



In a typical experiment, a pinhole (of radius R less than 
the Thomas Fermi radius) is often introduced just after 
the atomic cloud to restrict delay time detection to only 
light passing through the most dense central axial col- 
umn. This can be included by the averaged delay time 



td(T) 



ttR 2 



2ndr 



dz—-r 

l v g (r) 



2L 

c 



(14) 



For small pinhole size and long optical path length L, 
we can ignore the spatial dependence of L due to shape 
of the cloud |fl3|| . The group speed is then calculated 
according to v g — L(T)/td(T). Typical reports of our 



calculation are reported in Fig. [I]. Below the BEC tran- 
sition temperature we see that the group velocity in a 
Bose gas drops sharply. Deviations of Fermi gas from 
the classical behavior, on the other hand, appear below 
half of T F . For the parameters used N = 3.8 x 10 6 JO] 
and other parameters as given in Fig. [I] caption, Tp [31] 
is approximately twice the T c . In such degenerate low 
temperature regimes, different quantum statistics of the 
gases result in qualitatively and quantitatively different 
behaviors of their slowed group velocities. For this set of 
experimental parameters, we found local field correction 
causes only ~ 3% reduction of group velocity, invisible in 
the figure. 

The absorption caused by the imaginary part of the 
susceptibility (x") is best described by the the transmis- 
sion coefficient 
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which in turn gives the transmission T(A,T) = e Qr . A 
good estimate is provided by T = e~ 2 ^ oX ™ L ( T \ with Xm 
the imaginary part at the peak density (at r — 0) . Our 
numerical results are presented in Fig. ^ for T ~ T c /2 s» 
Tp/4, in the quantum degenerate regime. We notice 
the transmission for a Bose gas is always lowest, while 
for the corresponding Fermi is always the highest. This 
can be simply explained in terms if their peak densities. 
For near-resonant pulse propagation with |A| < 37 Bose 
gas becomes essentially opaque, while significant trans- 
mission through a Fermi gas is still possible. This high 
transmission feature can be considered as an advantage 
of a Fermi gas, over the Bose gas regarding applications 
in quantum memories, and optical data storage. 

It is interesting to discuss briefly results for the zero 
temperature limit, in which analytical solutions exist for 
all density distribution functions, given respectively by 

PF(T = 0) = -^(^-r 2 -eV)^, 

7T rip 

^ r = »> = S (B| -*- S - eV) ' 

where R B = (15N ea sc /a ho ) 1/5 a r and R F = (48iVe) 1/6 a r 
are the Thomas-Fermi radius for bosonic and fermionic 
clouds and a r = y / h/(Muj r ). 

After some tedious calculations, we obtain 
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for bose gas and 
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for fermi gas respectively at T = 0. 



Valuable insight can be gained by examining the ex- 
ternal parameter dependence of the above two formu- 
lae. Most notably v g ~ iV~ 2 / 5 for a Bose gas while 
v g ~ N^ 1 / 2 for a Fermi gas. Larger delay times are 
obtained for R <C i?s than for i? ~ i?p , Rb, and de- 
pends on the statistical nature of the atomic medium. In 
a Fermi gas, v g is three times larger for R ~ Rp than for 
i? -C Rf- While in the case of Bose gas, v g is 2.5 times 
larger for R ~ i?g than for i? <C This dependence 
may be used to calibrate and determine experimental ob- 
servables such as the cloud size, temperature, and the 
Thomas- Fermi radius. 

In conclusion, we have examined the temperature de- 
pendent behavior of slow light transmission through an 
ultra-cold quantum gas of two level atoms. We have in- 
cluded local field corrections, and compared the results 
for atoms with different quantum statistics. Analytical 
results are obtained below the condensation and Fermi 
temperatures. It is interesting to point out that as dis- 
played in Fig. |], the logarithmic of group velocity for a 
Fermi gas has a simple temperature dependence; it is al- 
most linear. This feature may be used for precise temper- 
ature measurements in regimes below the Fermi tempera- 
ture when precise calorimetry is difficult to achieve. The 
absolute scale can be obtained by a calibration against 
other independent measurements at higher temperatures. 
Furthermore, as emphasized in the theoretical studies 
in Ref. ]i~3[ ], a sharp discontinuity shows up for a Bose 
gas exactly at the condensation phase transition point. 
This is clearly displayed in Fig. |l|. For a three-level 
system, theoretical results for T < Tc are less reliable 
because of complications from both strong mean-field 
mutual couplings and multiple scattering effects. Con- 
sequently, agreement with experiment in this regime is 
poor pa]. In the two-level case studied here, our theory 
includes the self-mean-field interaction, and one may ex- 
pect its validity even in this ultra-low temperature regime 
as long as the incident probe field is weak. Unfortunately, 
there is no current experimental data to compare with. 
Should this prove to be true, this dramatic feature at the 
transition temperature can be used for accurate determi- 
nation of Tc- Finally we note that although these results 
are similar to those in three level EIT systems. The two 
level model studied here may present some practical ex- 
perimental advantages. 
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FIG. 1. Temperature dependent group velocity in Fermi 
(dash-dotted line), Bose (solid line), and Boltzman (dashed 
line) for two level trapped quantum gases. We used Na as 
an example with the following parameters: N — 3.8 x 10 6 , 
R = 7.5 (jum), w r = (2vr) 69 (Hz), e = 1/3, 7 = (2tt) 10.03 
(Hz), a sc = 2.75 (nm), w = (2vr) 5.1 x 10 14 (Hz) as in Ref. 
[1]. For these parameters Rb = 17.76 (/im), Rf = 50.04 
(jim), a r = 2.52 (fim), and a^ = 3.03 (/an). We chose the 
detuning A = IO7 to assure sufficiently transmission. 




FIG. 2. The same as in Fig. 1, except for the detun- 
ing dependence of the transmission coefficient. All three 
quantum gases are compared for the same temperature 
T = Tc/2 w Tf/4. 



